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Abstract 

Non-vascular plants (lichens and mosses) are significant components of tundra 
landscapes and may respond to climate change differently from vascular plants affecting 
ecosystem carbon balance. Remote sensing provides critical tools for monitoring plant 
cover types, as optical signals provide a way to scale from plot measurements to regional 
estimates of biophysical properties, for which spatial-temporal patterns may be analyzed. 
Gas exchange measurements were collected for pure patches of key vegetation functional 
types (lichens, mosses, and vascular plants) in sedge tundra at Barrow AK. These 
functional types were found to have three significantly different values of light use 
efficiency (LUE) with values of 0.013±0.001, 0.0018±0.0002, and 0.0012±0.0001 mol C 
mol 1 absorbed quanta for vascular plants, mosses and lichens, respectively. Discriminant 
analysis of the spectra reflectance of these patches identified five spectral bands that 
separated each of these vegetation functional types as well as nongreen material (bare 
soil, standing water, and dead leaves). These results were tested along a 100 m transect 
where midsummer spectral reflectance and vegetation coverage were measured at one 
meter intervals. 

Along the transect, area-averaged canopy LUE estimated from coverage fractions 
of the three functional types varied widely, even over short distances. The patch-level 
statistical discriminant functions applied to in situ hyperspectral reflectance data collected 
along the transect successfully unmixed cover fractions of the vegetation functional 
types. The unmixing functions, developed from the transect data, were applied to 30 m 
spatial resolution Earth Observing- 1 Hyperion imaging spectrometer data to examine 
variability in distribution of the vegetation functional types for an area near Barrow, AK. 
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1 Spatial variability of LUE was derived from the observed functional type distributions. 

2 Across this landscape, a fivefold variation in tundra LUE was observed. LUE calculated 

3 from the functional type cover fractions was also correlated to a spectral vegetation index 

4 developed to detect vegetation chlorophyll content. The concurrence of these alternate 

5 methods suggest that hyperspectral remote sensing can distinguish functionally distinct 

6 vegetation types and can be used to develop regional estimates of photo syntehtic LUE in 

7 tundra landscapes. 
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Introduction 

High northern latitudes are undergoing dramatic changes in climate. Warming 
trends have been observed in northern regions and Global Circulation Model predictions 
indicate arctic and boreal regions are likely to warm by several degrees over the next 
century, a greater rate than other regions on the globe [1, 2]. In high latitudes seasonal 
shifts in surface properties involve freezing and thawing of water, so relatively small 
temperature changes around the freezing point can result in significant environmental 
shifts. Along with temperature changes, significant changes in precipitation and 
evapotranspiration are also predicted for these regions [1,2]. 

Tundra vegetation is expected to be particularly responsive to climate change. 
Increased warming during the growing season is likely to alter production by lengthening 
the growing season and increasing metabolic activity [3, 4, 5]. Changes in seasonal air 
temperature may increase soil temperatures and the active layer depth to permafrost. 
This, in turn, affects soil microbial activity, nutrient cycles and soil moisture [6, 7], which 
alter existing plant growth and competitive species interactions resulting in dramatic 
changes in vegetation composition. Further, tundra ecosystems are often moisture 
limited, so climate induced changes in precipitation patterns and surface hydrology will 
also act to alter vegetation growth patterns [8,9]. 

A unique characteristic of the tundra ecosystem is the relative dominance of non- 
vascular plants in the landscape, such as mosses and lichens. In some cases, mosses 
dominate local habitats in the tundra because of their tolerance to extreme cold, 
dessication, water-logging, and low light [10]. Lichens can survive extreme climate 
conditions such as drought and persistent freezing as well [11, 12, 13]. Tundra warming 
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manipulation experiments have shown a decline in lichens with an increase in vascular 
plants suggesting changing proportions of tundra plant types as climate changes [14]. 

Non-vascular plants can provide a significant fraction of tundra carbon uptake and 
should be explicitly included in descriptions of tundra carbon fluxes. Photosynthesis of 
tundra vascular plant canopies has been described with strong relationships found 
between leaf area or the fraction of absorbed photosynthetically active radiation and gross 
primary productivity [15, 16, 3]. However, the net primary productivity (NPP) of mosses 
have been shown to represent about 25-30% of total above-ground NPP in several tundra 
sites including; tussock tundra [9], coastal tundra [17] and tundra heath [18]. In other 
cases mosses were found to dominate tundra photosynthetic C02 fluxes [19, 20]. 

There are further important physiological differences between vascular and non- 
vascular plants that affect tundra water and energy balance. Mosses and lichens do not 
have roots, so cannot access soil moisture below the surface. Instead they depend on 
water from atmospheric humidity, ground water, or precipitation. They also do not 
possess stomata and therefore lose water readily to the atmosphere and many species are 
adapted to survive long periods of desiccation. [21]. They may also act as an insulating 
layer to heat transport between the soil and the atmosphere [22]. Thus, in areas where 
these non-vascular plants make up a significant portion of the vegetation cover the 
response of ecosystem carbon and energy balance to environmental changes is expected 
to differ from that of vascular plant-dominated vegetation and must be accounted for in 
tundra ecosystem modeling. 

The difficulties in working in remote tundra locations have always incurred large 
logistical costs, therefore observations of ecosystem change using remote sensing provide 
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an expedient and economical method of collecting repeatable and consistent 
measurements over large areas. Importantly, most remote sensing theory of vegetation 
has been derived from studies focused on temperate crops and forests and thus the 
interpretation of conventional remote sensing tools is often ambiguous for northern 
landscapes. There are several unique features of the tundra that affect the interpretation 
of remotely sensed data [23]. Non-vascular plants can represent significant fractions of 
the tundra landscape cover, and these non-vascular plants have different spectral 
characteristics as well as different physiological responses from vascular plants or from 
bare soils [24, 25, 26]. For example, remote sensing studies have been able to detect 
changes in lichen-dominated areas due to responses to both short-term temperature 
anomalies as well as long-term temperature trends [27, 28]. 

This study investigates the importance of non-vascular plants in the tundra 
landscape and examines their effects on tundra carbon uptake, using spectral reflectance 
to distinguish these different functional types. We explore the concept of “optical types” 
[29] to scale from ground-based in situ measurements to landscapes using satellite 
observations to examine vegetation functional type patterns across the tundra. 

Materials and Methods 

Study Area 

The study examines an area mainly located north and west of the town of Barrow, 
AK, including most of the Barrow Environmental Observatory. Specific study plots are 
located at 71° 19' 19.1" N 156° 36' 15.9" W, approximately seven kilometers east of 
Barrow. The area is classified as moist acidic coastal tundra [30]. Barrow has a mean 
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annual temperature of -12.0°C, with a minimum mean monthly temperature of -26.6°C 
in February and maximum mean monthly temperature of 4.7°C in July. More than half of 
the 106 mm of annual precipitation falls as rain during a 3-month period from July 
through September, with the ground being snow covered from September through May to 
mid- June [31]. 

The tundra community at Barrow, AK consists of an overstory dominated by 
vascular plants (mostly graminoids) and an underlying mat of mosses [32]. Vascular 
plants at Barrow include dwarf shrubs, forbs, and perennial herbs but the most dominant 
growth form is graminoid [30]. Common vascular species are: Carex aquatilis, Dupontia 
fischerii, Eriophorum angustifolium, E. scheuzerii, Luzula confusa, Petasites frigidus, 
Potentilla hyparctica, and Salix rotundifolia. Mosses are a significant portion of the 
community and in some areas may account for more biomass than vascular plants [33] . 
Moss species include: Dicranum elongatum, D. undulatum, Drepanocladus revolvens, 
Polytrichum (juniperum), and Sphagnum spp. Lichens are scattered throughout the drier 
areas but are not as abundant or productive as mosses and graminoids [34] and include 
the species: Alectoria nigricans, Cetraria cucullata, C. nivalis, and Dactylina arctica. 

Measurements 

This project included: 1) the collection of in situ measurements of plots for the 
three vegetation functional types, including photosynthesis and spectral reflectance; 2) 
measurements of ground cover amount and spectral reflectance at every meter along a 
100 m transect; 3) development of relationships between spectral reflectance and both 
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photosynthesis and cover amount; 4) estimation of LUE across a tundra transect; and 5) 
application of these findings to the spectral information collected from Flyperion. 

Vegetation measurements were designed to relate photo synthetic rate to spectral 
reflectance. Tundra vegetation were grouped into three functional groups; vascular 
plants, mosses, and lichens, with at least four replicates per group measured on any one 
day. Each replicate provided a different species from each group in an attempt to 
measure the variability within a group. Field sampling for this study was conducted on 
July 20 and August 5, 2001, near the peak of the growing season in this area. 

Spectral reflectance measurements were collected at 51 plots, each about 15 cm in 
diameter, representing nearly pure areas of representative tundra species. Photo synthetic 
gas exchange measurements were collected at 17 of the plots. All of the study plots were 
in very close proximity (<15 m) to a 100 m linear transect that was also optically sampled 
at frequent intervals throughout the summer [35] . 

C0 2 gas exchange was measured at 17 plots where a 14.6 cm diameter collar of 
acrylic plastic was embedded in the ground. When measurements were made, the plot 
was enclosed in a clear acrylic chamber constructed of 0.3175 cm thick acrylic tube glued 
to a 0.635 cm thick acrylic top and a closed gas exchange system was connected to the 
chamber (LI-6200, LiCOR, Lincoln, NE). Each individual plot was measured under 
ambient environmental conditions for net photosynthesis and then the measurement was 
repeated with the chamber covered with a black cloth to block out light to provide an 
estimate of dark respiration. Gross photosynthesis was calculated as the difference 
between the net photosynthesis, measured with the clear chamber, and dark respiration, 
measured with the chamber covered. A quantum sensor mounted on the chamber top 
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measured incident photosynthetically active radiation (PAR). The chamber had a PAR 
transmittance of 95% [36]. Photo synthetic light-response curves were developed using 
household window screen material as neutral density filters to control incident PAR, 
allowing photosynthesis measurements for each functional type to be rapidly collected 
under a range of light levels. 

Light use efficiency (LUE) was determined as the ratio of absorbed light and 
gross photosynthesis: 



( 1 ) 


2 1 

where e is LUE, G is gross photosynthesis with units of ja mol CO 2 m " s , Qi is incident 

PAR with units of p mol photons m “ s (i.e. photosynthetic photon flux density, or 
PPFD) and /par is the fraction of the PAR absorbed by the vegetation (unitless). 
Absorbed PAR (APAR) is the product of /par and incident PAR. 

/par is not well defined and difficult to determine for vegetation with low growth 
forms such as lichens and mosses. We assumed that no PAR was transmitted through 
moss and lichen mats, so for these vegetation types /par is 1-PAR albedo. PAR albedo is 
estimated by integrating the measured spectral reflectance (see below) between 400 and 
700 nm. For the vascular plant plots /par is determined using the linear relationship 
between Normalized Difference Vegetation Index (NDVI) and /par described in 
Huemmrich et al. [36]. 

Optical Sampling 
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The optical field sampling consisted of measurements of the small single species 
plots, as well as measurements collected at every meter along a 100 m transect to observe 
areas of mixed coverage [35, 37, 38]. Optical measurements were collected using two 
portable field spectrometers (UniSpec, PP Systems, Haverhill, MA). These 
spectrometers measured reflectance between 310 and 1130 nm sampling at approximately 
3 nm intervals. 

For plot measurements, the spectrometer was fitted with a glass fiber optic cable 
connected to a stainless steel ferrule, which provided a 20° field of view. Reflected 
irradiance measurements were collected with the end of the fiber optic cable held 
vertically between 30 and 40 cm above the surface to view an area of less than 15 cm 
diameter to match the area of gas exchange measurements. Within seconds of collecting 
a surface measurement, a reference measurement was made of a calibration panel 
(Spectralon, Labsphere, Inc. North Sutton, NH). The speed of this procedure allowed the 
acquisition of data even under cloudy conditions. Three reflectance spectra were 
collected for each sampled plot, and were averaged to produce one spectrum per plot. In 
addition to the vegetation plots, spectral measurements for nongreen materials including 
bare soil, open water, and areas of standing dead vascular plants, were collected soon 
after snowmelt and before the start of green-up. For the vegetation plots, ground-based 
spectral reflectance measurements were collected near the height of peak productivity 
period in early August (August 5 and 8, 2001) and reflectance for nongreen plots were 
collected throughout the summer (June 20 and 30, July 15, August 8 and 22, and 
September 3, 2001). 
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To characterize spectral reflectance patterns at multiple spatial scales, spectral 
measurements were made at 1 m intervals along a 100 m linear track system [35, 38]. 
The track was supported less than a meter above the tundra surface by tripods with a tram 
cart riding on the track carrying a dual channel spectroradiometer (Unispec DC, PP 
Systems). From the spectrometer were two fiber-optic cables mounted on a mast and 
boom attached to the cart, one attached to a diffuser head viewing upward, the other on 
the boom pointing vertically downward to view an area south of the track with a field-of- 
view of approximately 1 m in diameter. With each measurement of reflected radiance, a 
measurement of incident irradiance was also acquired at the same time for a 
determination of surface reflectance. Due to the coincident measurements of both 
incoming and reflected radiation this system was able to correct for effects of varying 
cloud cover. The tram system collected spectral measurements that were repeatable 
throughout the growing season while minimizing disturbance to the tundra surface [35, 
38]. In this study we used measurements from the peak of the growing season period 
collected on August 10, 2001. 

Processing the ground-measured spectral reflectance data included interpolating 
the spectral bands to 1 nm intervals. Due to instrument noise at the ends of the 
spectrometer range, the usable spectral range was between 400 and 1000 nm. The 
processing used freely available software (Multispec v.4.0, 
http://specnet.info/specnet toolkit.htm) . The spectral reflectance data were then 
convolved to match the approximately 11 nm band passes of the Hyperion instrument 
[39] in subsequent analyses. 
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Imaging spectrometer data were acquired on July 20, 2009 (day 201) by the 
Hyperion instrument on the Earth Observing- 1. Hyperion provides data in 220 spectral 
bands from 0.4 to 2.5 pi m at 30-meter resolution and a 7.5 km wide swath [39]. The 
satellite data were atmospherically corrected to surface reflectance using the Atmosphere 
Removal algorithm (ATREM) [40, 41]. To match the spectral range of the ground-based 
spectral data only 55 of the continuous Hyperion bands between 437 and 993 nm were 
used. 

Discriminant analysis was used to separate the functional groups based on optical 
properties, determining the important spectral bands and simplifying the number of 
independent variables [42]. The plot data were divided into four groups: vascular plants, 
mosses, lichens, and nongreen materials (bare soil, dead vegetation, and standing water). 
Inputs to the discriminant analysis using Systat (Version 13, Systat Software, Inc., 
Chicago, IL) were plot reflectance spectra convolved to the Hyperion bands. The 
analysis used a forward stepwise method (Wilks’ Lamda) to test for significant mean 
group differences, using a probability of F< 0.05 to include bands and F> 0.10 to remove 
the band. 

At each meter along the 100 m transect a visual estimate of species abundance 
was made over the period August 8-11, 2001 using a 1 m by 1 m quadrat. The species 
data were converted into percent area of each of the vegetation functional types and 
nongreen material types. We estimate errors in the estimation of coverage to be +/- 15% 
cover. Because of the three dimensional structure of the vegetation the sum of the 
coverage was more than 100%. For this analysis, the total measured cover values were 
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normalized so total cover equals 100% as is generally done in linear unmixing 
approaches. 

LUE was determined for areas with mixtures of different functional types as the 
weighted sum of the LUE from the pure plots, where the weight coefficients were the 
normalized fractional cover values. This calculation assumes that the incident PAR is 
absorbed by the landscape components proportional to their normalized cover fraction 
and was applied to determine LUE for both the transect data and the Hyperion satellite 
imagery. 

Chlorophyll concentrations are a key physiological factor related to maximum 
photosynthetic rates. We estimated chlorophyll concentrations using Gitelson’s three- 
band model: 


where the spectral index C, is proportional to the chlorophyll concentration, and p ? is the 
reflectance for a given wavelength band, 7. Wavelengths chosen are maximally sensitive 
to absorption by chlorophyll and other pigments (7 1 =549 nm), are influenced primarily 
by non-chlorophyll pigments (72=793 nm), and where reflectance is controlled by leaf 
scattering (73=793 nm) [43, 44, 45]. 

The cover estimates of the functional types retrieved from the Hyperion imagery 
were used to derive landscape LUE patterns. The LUE per pixel was calculated as 
described above and compared with chlorophyll concentration estimates obtained using 
Equation 2. 



( 2 ) 


Results 
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The physiological measurements identify differences relating to photo synthetic 
carbon exchange among the three vegetation functional types. From the plot gas 
exchange data, LUE for each functional type was calculated as the slope of a linear 
regression between gross photosynthesis and APAR, forced through the origin (Figure 1, 
Table 1). Lichens, mosses, and vascular plants all had different LUE values. Vascular 
plants had significantly higher net photosynthetic rates and LUE than the other two 
functional types. The vascular plant LUE from this study (0.0134 mol C mol 1 quanta) 
was close to that previously determined for a vascular plant overstory near Barrow 
(LUE=0.0126 mol C mol 1 quanta [36]) and within the 95% confidence interval of the 
slope from the regression in this study (i.e. the vascular plant LUE). This was noteworthy 
as the LUE in this study was based only on measurements collected during the middle of 
the growing season, while the Huemmrich et al. [36] LUE was derived from 
measurements collected throughout the growing season, suggesting the LUE for vascular 
plants was relatively stable over much of the season, as well as having similar LUE 
values for vascular plants in mixtures with mosses or growing alone. 

The LUE for mosses and lichens were much lower than for vascular plants. There 
was a significant difference, however, between these two functional types (Table 1) with 
mosses having a LUE almost 50% greater than lichens (a similar analysis for LUE of the 
moss data was presented in [36]). 

The in situ spectral reflectance measurements of the plots were divided up into 
four groups: vascular plants (n=22), mosses (n=ll), lichens (n=18), and nongreen 
material (n=19) (Ligure 2). Generally, the vascular plant and moss spectra were similar, 
and display typical green plant spectral reflectance patterns, including more variation in 


13 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


Functional Types and Unmixing 


3/18/13 


the near infrared than in visible wavelengths. In comparison, lichen spectra had much 
higher visible reflectance and were more variable in all wavelengths. Spectral patterns for 
the nongreen materials were clearly different for dead vegetation, bare soil, and water. 
Generally the spectral variability for the entire nongreen group increases as a function of 
wavelength. 

Discriminate analysis provided an objective tool to evaluate the ability of optical 
sampling to distinguish functional groups and to identify a subset of the spectral bands to 
do this separation. The stepwise method for discriminate analysis reduced the number of 
wavebands from 55 to five, yet retained separation of functional types. The five chosen 
bands were located in key locations of typical green plant spectra (Figure 2); at the blue 
(488 nm) and red (671 nm) chlorophyll absorption wells, near the inflection point of the 
red edge (712 nm), at the shoulder of the red edge (763 nm), and in the near infrared 
“plateau” (834 nm). Mulhern [46] identified the blue band as an important spectral region 
for separating lichens from soils. The discriminate analysis functions (Figure 3) correctly 
predicted 83% of group membership overall, with 86% correct for vascular plants, 91% 
for mosses, 95% for the nongreen materials, and 61% for lichens. Lichen errors were 
mainly due to confusion with the nongreen materials (Figure 3). 

Even at the scale of one square meter, the cover estimates along the tram transect 
show that all square meter plots consisted of mixtures of multiple cover types (Figure 4). 
The observed cover estimates included little bare ground with maximum coverage of 
20%, while standing dead vegetation ranged from 0-70%, and water from 0-100%. Moss 
coverage ranged from 0-100%, but when normalized to a 100% maximum total coverage 
maximum moss coverage was 40% +/- 7%, based on the propagation of an assumed 15% 
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error in the original observations. Lichen coverage ranged from 0-90% with 35% +/- 7% 
maximum normalized lichen coverage. Vascular plant coverage ranged between 40- 
150%, with normalized coverage of 18-88% +/- 9%. Vascular plant coverage was highest 
in the locally low areas, particularly the areas with standing water, while moss coverage 
was highest in locally low areas without standing water and on the edges of the wet areas. 
Lichen coverage was highest in the locally high areas where bare soil patches also 
occurred. These spatial patterns indicate the role of microtopography on the distribution 
of the functional types. 

The transect coverage observations indicate how heterogeneous the tundra is, as 
areas over a few square centimeters in size generally consist of mixtures of multiple 
types. Therefore, to estimate LUE at that scale of a 30 m Hyperion pixel, the fractions of 
each functional type must be determined. We used the optical measurements to scale 
from plot-level LUE to 30 m areas. Since the discriminant functions optimized separation 
of the different functional types with a small number of spectral bands, we used the 
distance to endmembers in the discriminate analysis function space to estimate cover 
fractions in mixed pixels. Assuming that the range of in situ plot reflectances provide a 
reasonable description of the variability of plot type, the endmembers for the moss, 
lichen, and nongreen cover types were assigned values of the average of the plot 
reflectances. However, variability in the vascular plants reflectance is driven by 
differences in green leaf area index. To account for that variability the vascular plant 
endmember was chosen to be the “greenest” plot spectra (i.e. spectra with the highest 
NDVI, which was also the greatest distance from the average of all the plot values in the 
discriminate analysis function space). 
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Endmember reflectance for each functional type and reflectances at each meter 
along the transect were transformed using the discriminate analysis functions. The 
distances in the discriminate analysis space between the coverage estimates for each 
square meter and the endmembers were calculated, these were related to the observed 
coverage to create transfer functions. The statistical associations between distance from 
endmembers and observed coverage were low due to large observation errors relative to 
the range of values, with R 2 values of 0.39, 0.13, and 0.19 for vascular plants, mosses, 
and lichens respectively. Nevertheless, the cover fractions determined from the in situ 
reflectance were found to be in close agreement with sampled cover fractions (Figure 6), 
with root mean square error (RMSE) of 11% cover for vascular plants, 8% for moss, 9% 
for lichen, and an overall RMSE of 9% cover. 

LUE was calculated for each square meter block along the transect using the 
observed cover fraction and coverage estimated from the reflectance spectra (Figures 5 
and 7). Along the transect LUE was quite variable, with maximum and minimum values 
occurring within 6 m of each other (Figure 7). Variability in LUE was affected by 
microtopography, with higher values of LUE tending to occur in locally low, wet areas. 
This variability over the 100 m distance produced a fourfold difference in LUE with 
values from 0.003 to 0.012 mol C mol 1 quanta for the in situ observations. The average 
errors for estimating LUE were similar for the two methods: 0.0014 mol C mol 1 quanta 
for the method using observed coverage and 0.0015 mol C mol 1 quanta for the LUE 
based on the remote sensing approach. The two different LUE calculations have a 
correlation of 0.62 with a RMSE of 0.0013 mol C mol 1 quanta. 
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The equations relating spectral reflectance and functional type coverage derived 
from the transect data were the applied to the Hyperion imagery to create a continuous 
fields description of the coverage of the functional types for the tundra around Barrow 
(Figure 8b). Over this area the distribution of vascular plant coverage per pixel had a 
maximum value of 74% with a distribution peak of 32%, for mosses the maximum 
coverage was 31% with the distribution peak at 22%, and for lichens the maximum 
coverage was 22% with a peak of 10%. 

Without ground observations of vegetation type coverage, an examination of an 
area known to be free of vegetation was used to make an error estimate. The old Naval 
Arctic Research Laboratory runway was chosen as for this test. For runway pixels 
vascular plant coverage was as high as 9%, moss 14%, and lichen 13%, compared to 
expected values of zero percent coverage. These values were close to the RMSE of 
coverage from the tram analysis comparing the spectral unmixing with ground 
observations, suggesting the Hyperion unmixing captured the actual functional type 
distribution near this level of error. 

There are clearly observed spatial patterns in the cover fractions in this landscape. 
High vascular plant coverage occurred on the margins of drained lakes. Higher lichen 
coverage tended to be found in drier upland regions. This pattern of lichen occurrence 
was also observed in the transect measurements (figure 4). Moss coverage was more 
widely distributed, although high moss coverage in the middle of the drained lakes may 
be erroneous due to confusion between mosses and mixtures of water and vascular plants. 
Field observations confirm that moss often occurs as a low, “background” layer beneath 
the vascular canopy of wet tundra in this region [36] . 
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The Hyperion functional type coverage estimates were used to calculate landscape 
patterns of LUE (Figure 8c). For the vegetated parts of the scene, this Hyperion-based 
FUE ranged from 0.0021-0.0102 mol C mol 1 quanta, a fivefold difference between 
maximum and minimum FUE. The distribution peak of FUE for this area and the overall 
average were both 0.0048 mol C mol 1 quanta. 

The Hyperion-based FUE was compared with the chlorophyll spectral index, Ci 
(Equation 2). Ci was well correlated with FUE over the study area (R 2 =0.69, SE=0.0006 
mol C mol 1 quanta) (Figure 9) suggesting vegetation chlorophyll concentration is a key 
determinant of FUE for this tundra ecosystem. This suggests that a more direct approach 
for deriving spatial patterns of FUE for tundra from imaging spectrometry would be 
based on a simple reflectance-based metric of chlorophyll concentration. 

Discussion and Conclusions 

This study uses the concept of optical types [29] to examine a key characteristic 
of ecosystem carbon exchange in the tundra. To have usable optical types one must be 
able to both identify significant functional differences in vegetation types and spectrally 
identify these different types. Our grouping of tundra vegetation (vascular plants, mosses, 
and lichens) meets these criteria. Among the different types, the plot measurements 
showed distinct differences in FUE, a key variable describing photosynthetic carbon 
uptake and discriminate analysis showed that the types could be separated based on their 
spectral reflectance. The ability to link FUE to reflectance characteristics provides an 
approach to scale from the ground measurements to distributions over a landscape. Using 
an unmixing approach we were able to create continuous fields of key tundra functional 
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types for this area, an improvement over simple classifications. Results from this study 
demonstrated how, even at a scale of a square meter, variable mixtures of functional 
types produced significantly different values for LUE. Unmixing using spectral 
reflectance provides a unique estimate of LUE for each pixel. Remote sensing provides a 
way to map large areas and make noninvasive repeat measurements to monitor ecosystem 
change. Remote sensing is a particularly important tool for observing tundra due to the 
difficulties and expense involved in traveling to and working in this region, rendering 
direct field sampling over large areas impractical. 

The linear unmixing approach used to calculate pixel-level LUE simplifies the 
tundra canopy as it does not take into account the three dimensional structure of a 
vascular plant overstory. Vascular plants can grow through mats of mosses and lichens 
and shade them, decreasing their photo synthetic production [36] . However, vascular plant 
leaf area index for this region is generally low, with observed midseason values less than 
2 [36] . We believe the overall effect of this overstory on the moss and lichen understory 
production is relatively small in this region due to both the low leaf area combined with 
generally diffuse light conditions. 

This study provided a snap-shot of this landscape at a particular point in time. 
Seasonal growth of vascular plants will change their relative proportions throughout the 
growing season and this affects carbon uptake and LUE patterns [38]. Early in the 
growing season, as vascular plants begin to grow, tundra productivity will be dominated 
by moss and lichen photosynthesis, with the vascular plant component increasing as the 
season progresses [36]. Thus, the spatial patterns of coverage and LUE are expected to 
change through the growing season. 
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The correlation between LUE calculated from coverage estimates and the 
chlorophyll index (Ci), even though they use different algorithms and spectral bands, 
demonstrates how multiple approaches can yield similar results. This convergence is 
most likely due to plant physiological constraints indicating chlorophyll concentration is 
a controlling factor in vegetation photo synthetic efficiency, and may provide a more 
direct approach to deriving LUE. 

This study found significant variability in functional type cover fractions, leading 
to variability in estimated LUE, at local scales (meters) in the field data collected along 
the transect and at landscape scales (km) in the satellite imagery. The spatial variability 
shown in LUE estimates is not accounted for in existing carbon flux models. Spatial 
distributions of vegetation functional types were strongly related to surface 
microtopography, from which we infer surface hydrology to be a significant controlling 
factor. These results suggest that climate change affecting surface hydrology [47] will 
also affect spatial patterns of vegetation distribution and ecosystem carbon exchange. 
Particularly in a scenario of modified surface hydrology, we can expect that climate 
change will affect relative coverage due to differing responses of each of the functional 
types to climate change, and that these changes in cover will be associated with altered 
carbon and energy balance. 

Since hydrology is an important factor determining the vegetation functional type 
coverage, shortwave infrared bands may provide additional information on vegetation 
water content, improving the retrieval of functional type coverage [48]. This approach 
was not tested in this study as the field measurements observed only visible and near 
infrared spectral bands. Recent observations combining NIR and blue bands to detect 
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standing water in coastal tundra [49] may also prove useful at detecting hydrological 
changes linked to cover type and LUE. 

The results of this study are indicative only of this part of the Alaskan coastal 
tundra and further field studies are required in other parts of the tundra biome to test and 
extend the results and to evaluate the ability to generalize the results. For example, shrubs 
are an important component of the tundra biome and shrub coverage has been shown to 
be increasing in the tundra [50]. Shrubs did not represent a significant cover type in our 
study area, however, so were not explicitly included in this study. Future work is required 
to examine shrub optical and physiological characteristics to see if they comprise an 
optically distinct functional type. 

The use of multiple narrow spectral bands for the umixing points to an application 
for the future NASA Hyperspectral Infrared Imager (HyspIRI) mission for mapping 
tundra vegetation distribution according to physiological function, hydrology, and 
microtopography over the entire biome. The concept of optical types (optically 
distinguishable functional types) may prove to be useful for this kind of functional 
mapping. Similarly, simple metrics of surface water cover (e.g. [49]) or pigment 
distribution (e.g. [45]) may provide insight into the proximal causes of functional change 
in tundra ecosystems. 

The close linking of optical properties and carbon flux measurements in this study 
has provided insights into the functioning of the tundra ecosystem by identifying fine- 
scale patterns of FUE that may be related to hydrology and microtopography. This study 
also illustrates new “scaleable” functional mapping techniques that can be extended 
through the use of satellite remote sensing to larger regions. This approach lends itself to 


21 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


Functional Types and Unmixing 


3/18/13 


future studies considering how species composition relates to changing ecosystem 
function and provides a framework for studying ecosystem change through remote 
sensing that considers shifts in hydrology, species composition, and their effects on 
carbon balance. 
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1 Tables 



Vascular 

Plants 

Moss 

Lichen 

LUE 

0.01337 

0.00184 

0.00124 

St. Err. of 
Slope 

0.00101 

0.00022 

0.00011 

Num. Obs. 

29 

25 

28 

R 2 

0.74 

0.56 

0.74 

P 

<0.01 

<0.01 

<0.01 


3 Table 1. Light Use Efficiency based on regressions between net photosynthesis and 

4 APAR from plot data. LUE is in units of mol C mol 1 absorbed quanta. 



Vascular 

Plants 

Moss 

Lichen 

R 2 

0.39 

0.13 

0.19 

P 

<0.01 

<0.01 

<0.01 

Est. Obs. Error 

9.5 

7.1 

7.0 

SE of Regression 

10.8 

8.4 

9.0 


6 Table 2. Description of relationships in spectral space between the statistical distance to 

7 endmembers and coverage estimate observations, acquired along the tram transect. The 

8 errors are expressed in percent coverage. For all cases the number of observations is 100. 


31 


Functional Types and Unmixing 


3/18/13 


1 Figures 


2 

3 


(A 

CM 

E 

CM 

O 

o 


o 

E 

3 

(A 

'(A 

d) 

SI 

+■> 

c 

> 

(A 

O 


(A 

(A 

O 

i_ 

0 



0 200 400 600 800 1000 1200 1400 1600 


Absorbed PAR (pmol nr 2 s 1 ) 


Figure 1. Absorbed PAR and gross photosynthesis for the three vegetation functional 


4 types from plot data collected July 20 and August 5, 2001. Light use efficiency for each 

5 functional type was calculated from these data as the slope forced through the origin. See 


6 Table 1 for description of regressions. 
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Figure 2. Endmember spectral reflectance for each cover type, vertical lines indicate 
wavelengths used in the discriminant analysis. 
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Figure 3. The distribution of data from plots using the first two functions derived from 
discriminate analysis for the three vegetation functional types and nongreen materials. 
The larger gray circles indicate values used as endmembers in the unmixing analysis. For 
the endmember points the symbol within the circle indicates the vegetation functional 
type of that endmember. 
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Figure 4. Observed functional type coverage at 1 m intervals along the transect, visually 
estimated over the period August 8-11, 2001 . Because the vegetation is 3-dimentional, 
total cover fractions were over 100%. In the data analysis the total cover fractions were 
normalized to equal 100%. 



8 Figure 5. LUE along the 100 m transect shown with the black line was estimated from the 

9 functional type coverage (shown in Figure 4) and light use efficiency (shown in Figure 

10 1), along with LUE estimated using coverage derived from spectral unmixing shown with 
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the gray line. Microtopography is shown as dashed line. Error in Observed LUE estimate 
is 0.00136 mol C mol 1 quanta, RMSE in LUE derived from observations compared to 
LUE from spectra is 0.00266 mol C mol 1 quanta. 
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Figure 6. Comparison between observed vegetation functional type coverage normalized 
to 100% for vascular plants, mosses, and lichens along the transect and cover fractions 
estimated from spectral reflectance. Line is 1 to 1 line. The absolute error in visual 
estimates is approximately 9%. The RMSE of spectral retrievals is 9% absolute. 
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2 Figure 7. LUE (mol C mol 1 quanta) along the transect estimated based on observed 

3 normalized cover fractions and cover fractions obtained from spectral unmixing. Line is 

4 1 to 1 line. Errors in LUE based on observed cover fractions are due to errors in LUE for 

5 the pure functional types and an assumed 15% error of observed raw coverage amounts. 

6 Errors in LUE based on spectrally derived cover fractions are due to errors in LUE for 

7 pure functional types and the RMSE of the spectral retrieval of cover fractions. 

8 


37 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 


Functional Types and Unmixing 


3/18/13 



Figure 8. Three versions of the EO-1 Hyperion image acquired on July 20, 2009 are 
shown. 

Left: 3-band (RGB=834, 671, and 549 nm) composite image of surface reflectance. The 
grid of light blue lines on the lower left side of the image is the city of Barrow. The 
straight blue line along the shore near the top of the image is the old airport runway, used 
in the error evaluation. The oblong features scattered around the region are drained 
thermokarst lakes and the dark red ones are now marshes. 

Middle: Three band RGB continuous fields of estimated coverage of vegetation 
functional types derived from spectral unmixing and scaled between 0 and 50% coverage. 
R = Vascular Plant Cover, G = Moss Cover, B = Lichen Cover 

Right image: Map of LUE spatial patterns (mol C mol 1 quanta xlOOO) based on coverage 
estimates. 
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1 

2 Figure 9. Scatter plot of chlorophyll index (x-axis) versus LUE (mol C mol 1 2 3 4 5 quanta 

3 xlOOO) based on normalized cover fractions (y-axis). Both variables were retrieved from 

4 the Hyperion image. 
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